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We present results for the leptonic decays of heavy-light mesons in the quenched approximation and a prelim- 
inary result for /d, with dynamical fermions for Nf = 2. We examine the systematic effects of extrapolating the 
decay constants from the charm to the bottom quark mass, and find this extrapolation to be under control. We 
present a preliminary result for /^^ and compare to /j,^ matched to have the same lattice spacing set by 
the Sommer scale, ro. We find no significant difference. 



1. INTRODUCTION 

The decay constants of the B and D mesons 
are important in phenomenology. The combina- 
tion /bV Bb, \/Bbs, and their ratio, play a 
crucial role in extracting CKM matrix elements 
from experimental data. The B-factory exper- 
iments BaBar and Belle will measure — 
mixing quantities, but measuring will be very 
difScult. 

On the lattice, /s is a mature calculation, and 
less extensive than Bb, so we can test the sys- 
tematics of a heavy-light lattice calculation. The 
"world average" for Jb 0] in the quenched ap- 
proximation is Jb = 175(20) MeV and with dy- 
namical fermions (DP), fg^^^ = 200(30) MeV 
(see also ||]). This suggests a strong dependence 
on sea quark effects. However, it is not clear how 
to disentangle different systematic effects. 

2. QUENCHED DECAY CONSTANTS 

These results have been partially reported in 
nil. We use the non-perturbatively (NP) 0{a) 
improved Wilson quark action at two values of 
the coupling, P = 6.2 and (3 = 6.0. The 
light quarks have masses around strange, and the 
heavy quarks around charm. In the NP improve- 
ment scheme, the axial current is given by 

- Za{1 + bAam^){A^ + acAd^P,) (1) 

We use non-perturbative determinations of the 
coefficients Za, bA and ca by Bhattacharya et 
al. 1^ . In 1^ we show that the value of the de- 
cay constant depends strongly on the value of the 
mixing coefficient ca- These coefficients are only 
defined up to 0{a), e.g. compare the ALPHA § 
and Bhattacharya et al. values of ca &t f3 ^ 6.0. 



^ALPHA ^ _o.083(5) and c^'^^" = -0.037(4)(7). 
However, observables computed with these coef- 
ficients have the same continuum limit, provided 
each set of coefficients is determined consistently. 

We use continuum Heavy Quark Symmetry 
(HQS) to motivate the form of the extrapolation 
to the b quark mass; 

$(M)=e(M)/A/M = 7(l + ^ + ^) (2) 

where M is the pseudoscalar meson mass, / the 
pseudoscalar decay constant defined by 

{Q\A^\Ps{p))^ifPf, (3) 

and 8 is the logarithmic correction from QCD. 

This extrapolation is shown in figure |l|. We 
have also computed the value of $ in the static 
limit, but have not included the improvement 
term. We note that ca is negative, which will re- 
duce the value of We do not interpolate from 
charm, but use the static point as a check in the 
extrapolation. 

Also shown in figure |l| is a change of normalisa- 
tion of the axial current labelled CB, as suggested 
by Bernard ||l|, 

«)cB =Za{1 + bAanig + 2acAnf'^ (4) 

where 7^ = do{Q\P\P) /mq{Q\Ao\P) This dif- 
fers from the normalisation in equation (^) at 
0(0^), but has a finite static limit. This is 
plotted against M2, the kinetic mass, Mj"^ = 
d^E/dpf\p_^ There is a small difference between 
our standard extrapolation and that labelled CB. 

To estimate the effect of mass dependent lattice 
artefacts directly from the data, we add terms to 
equation (^) that model the lattice artefacts, 

$'(Af, a) = $(Af ) + 7 (£(aM)2 + C{aM)^) (5) 
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Figure 1. HQS extrapolation of $ at /3 = 6.0 
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Figure 2. Quasi-Continuum extrapolation of <i> 



We then fit to both lattice spacings simultane- 
ously. Shown in figure || is the data at both values 
of (3. The fit parameters e and ( are small and 
consistent with zero. Line A(B) shows the func- 
tion $'(M, a) with a at /3 = 6.2(6.0). These show 
that lattice artefacts indeed do blow up if extrap- 
olated. However, the lines labelled "/3 = 6.2(6.0) 
only" and "quasi-continuum" (QC) ($'(M,0)), 
lie close together. Moreover /^^ = 186(12) MeV, 
is close to the values of fs at fixed lattice spac- 
ing, albeit with larger errors, which implies mass 
dependent lattice artefacts are not spoiling the 
extrapolation. 

We estimate the systematic uncertainties from 
the quantity used to fix a, (ro,mp), linear vs 
quadratic HQS extrapolation, quasi-continuum 
extrapolation and the difference between differ- 
ent lattice spacings. We observe good scaling be- 
tween different lattice spacings. The values for 
the quenched decay constants are as follows, with 
the scale set by . The errors are statistical and 
systematic; 



fB = l9m-lt MeV 
/b, =220(6)11 MeV 

^=1.13(l)ii 
fK=150{3)t'l MeV 



fD=206{4)tll MeV 
fDs^'229{3)tll MeV 



3. DECAY CONSTANTS WITH Nf = 2 

The UKQCD collaboration has a set of 
matched ensembles 0, where the lattice spacing 
(as measured by tq) is fixed to be the same, as 
P and rrig'''^ are varied, including a quenched set. 
These ensembles are NP C(a) improved. Mea- 
suring the decay constants on these ensembles 
will enable us to study the effect of unquench- 
ing whilst keeping everything else fixed. Further- 
more, we intend to study the systematics of differ- 
ent heavy quark methods by comparing NRQCD, 
FNAL formalism and charm plus extrapolation 
directly on the same data sets. In particular it is 
interesting to see if the ratio /^^~^//^^~° is the 
same for each method. 

We present a preliminary result for f]j_, at fixed 
lattice spacing, comparing the Nf = 2 value with 
the quenched value. The simulation details are 
shown in Table |l| 

Table 1 

Simulation details 



/ 



/3, K""" 

cfgs 
a^^ , volume 
valence light 
valence heavy 
CA BPT 
CA NP 



5.2, 0.1350 5.93, 
0.70(1) 1 
202 117 
0.1 fm, 16^ X 32 

3 ~ strange 

4 ^ charm 
-0.0164 -0.0132 
-0.04(2) -0.032(14) 



3 



The values of the improvement coefficients have 
been determined to one-loop in perturbation the- 
ory 1^ . At one- loop these are independent of sea 
quark effects. We evaluate these in "boosted" 
perturbation theory. The interpolation to charm 
is shown in figure We see no significant differ- 
ence. 
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Figure 3. Comparison between DF and quenched 
data 

There is a preliminary NP determination 
Table || shows the dependence of Jd, on the value 
of CA- The entry Ps shows the value of the pseu- 
doscalar density contribution. 

Table 2 

f]j^ in MeV and its dependence on ca- 



CA 


Nf ^2 


Nf = 


BPT 


253(9) 


253(6) 


NP 


237(9) 


241(6) 





264(9) 


262(6) 


Ps 


693(24) 


629(14) 



Other than a small difference in the improve- 
ment term there is no difference. One reason may 
be the sea quarks we have used are still rather 
heavy. The Sommer scale vq is a short distance 
quantity and the decay constant, being propor- 
tional to the wavefunction at the origin is also a 



short distance quantity. We have fixed a with a 
short distance quantity, thus the short distance 
physics will be similar for both quenched and un- 
quenched. It is not surprising that we observe 
no difference. The MILC collaboration ||l^ see a 
similar pattern. Little difference when the scale 
is set by ri (related to tq) and significant rrig'^'^ 
dependence when the scale is set by rup. The 
decay constants are dimensionfuU quantities, so 
they are are sensitive to the scale setting ambigu- 
ity in the quenched approximation. The change 
in fDs from unquenching depends very strongly 
on how we choose to set the scale. 
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